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We consider a model of multiband superconductivity, inspired by iron pnictides, in which three 
bands are connected via repulsive pair-scattering terms. Generically, three distinct superconducting 
states arise within such a model. Two of them are straightforward generalizations of the two- 
gap order parameter while the third one corresponds to a time-reversal symmetry breaking order 
parameter, altogether absent within the two-band model. Potential observation of such a genuinely 
frustrated state would be a particularly vivid manifestation of the repulsive interactions being at 
the root of iron-based high temperature superconductivity. We construct the phase diagram of this 
model and discuss its relevance to the iron pnictides family of high temperature superconductors. 
We also study the case of the Josephson coupling between a two-band s' or s± superconductor and 
a single-gap s-wave superconductor, and the associated phase diagram. 
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I. INTRODUCTION 

The discovery of a new high-temperature supercon- 
ducting family of iron-based materials^— and the subse- 
quent developments have brought the question of multi- 
band superconductivity to the forefront of the condensed 
matter research. First discussed fifty years ago, this 
problem had remained somewhat obscure until iron pnic- 
tides, with their multi-band Fermi surfaces, made its un- 
derstanding an intellectual imperative. Following the 
initial discovery, an order parameter based on a two- 
band model was proposed as a likely possibility^. Soon 
thereafter, this so-called extended s-wave (or s± or s') 
superconducting state has been shown to be favored 
by the multiband electron dynamics of iron pnictides, 
both within a random phase approximation (RPA) type 
picture^— and in various renormalization group based 
approaches^— - as well as arising from a strongly corre- 
lated local limi t 17 i 18 - and is currently viewed as the most 
plausible superconducting state for these compounds. 

The first theoretical studies of a multi-band 
superconductivit y 1 ^ 20 were a straightforward gen- 
eralization of the BCS theory, with gap equations for 
several bands and attractive interactions. The most 
interesting result was that the (two) superconducting 
gaps A'(T) do not necessarily follow the single-gap 
BCS temperature dependence. Soon, however, it was 
realized 21 that the two-band model brings something 
conceptually new - superconductivity can be enhanced 
even by purely repulsive interband interaction. This 
requires a relative minus sign between the gaps on 
different portions of the multiply connected Fermi 
surface, while otherwise retaining an overall s-wave 
symmetry. In this way, it was argued in Ref. l2ll . the 
electron-phonon superconductivity in transition metals 
could receive an additional boost from the Coulomb 
repulsion driven resonant pair scattering between the 
broad s or p bands and narrow d bands at the Fermi 
level. 

The above conceptual novelty, however, extends much 



deeper than anticipated in Ref. |21| . The purely elec- 
tronic interactions could, in principle, produce supercon- 
ductivity even in the absence of any phonon-mediated at- 
traction. The superconductivity in this case would arise 
solely through the resonant pair scattering between the 
two bands and both, or more as the case may be, of 
these bands could be narrow d or even f bands. This 
promising mechanism for achieving high-temperature su- 
perconductivity - using purely electron-electron interac- 
tions with cut-off of order of Fermi energy instead of De- 
bye frequency, and thus potentially much higher transi- 
tion temperature T c - remained, however, largely ignored 
for the next fifty years. The reason is, basically, that the 
conditions in real materials are less than favorable. For 
s'-state to be operational the superconductivity-driving 
interband pair scattering has to be stronger than the 
superconductivity-suppressing intraband repulsion (most 
commonly they both come from the screened Coulomb 
interaction in metals). This is unlikely for at least one 
reason: the interband interaction usually involves higher 
momentum transfer, bands typically being well-separated 
in the fe-space, and is generically smaller. Thus, the sign- 
changing order parameter was considered unrealistic. 

This perception changed last year, with the advent of 
iron pnictides. At least for those members of this fam- 
ily that exhibit the highest T c s, a nodeless multi-gap 
order parameter - with some ARPES experiments see- 
ing as many as four different gaps^ 2 - - appears firmly es- 
tablished. The conventional electron-phonon interaction 
seems too weak to explain T c as high as 57 K (although 
some highly unconventional strong electron-phonon cou- 
pling still remains a remote possibility). This state of 
affairs makes the purely repulsive electronic interaction 
as the source of superconductivity and the s'-wave state 
in particular very appealing, even though it is still not en- 
tirely clear how the generic repulsion problem, described 
in the previous paragraph, can be overcome. There are 
some very promising studies in this direction, based on 
the renormalization group arguments and the peculiar 



band structure of these compounds 
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suggesting a 
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plausible route to this superconducting state. 



Many theoretical studies so far have used some varia- 
tion of the two-band model. A number of useful results 
were derived and valuable insight was gained within this 
simplified picture 23 . The real materials, however, are 
more complex, and some tight-binding representations of 
iron-pnictides 1 ^ indicate that typically three bands - one 
electron and two hole-like - are those most strongly cou- 
pled in the pair-scattering channel (see also Ref. wM . All 
this adds some urgency to the study of multiband super- 
conductivity with three or four bands. In this paper, we 
concentrate on a rather generic three-band model with 
repulsive interactions. The main question we are inter- 
ested in is " Is something conceptually new emerging from 
this increase of the number of the bands?" . The answer is 
"Yes," despite the fact that the gap equations themselves 
have the appearance of straightforward generalization of 
the two-band case. The reason for this is the frustration 
which the additional band introduces into the problem. 



To develop some intuitive understanding of the model 
let us start with an effective two-band situation. We 
ignore the intra-band interaction and consider only iden- 
tical bands. If the coupling with the third band is neg- 
ligible, there are two gaps A 1 = —A 2 and the overall 
magnitude is determined by the standard BCS relation. 
If we now introduce coupling to the additional band there 
are several possibilities. The system can stay in a two- 
gap state - now there are three such states - and keep 
the remaining band (nearly) gapless. In that sense, the 
interactions between the bands are frustrated, i.e. with 
such superconducting order one of the bands will not 
achieve what would otherwise be its natural gapped state. 
As suggested previously 2 ^., there is also a possibility for 
a new superconducting order parameter which compro- 
mises between the different frustrated two-gap order pa- 
rameters. We show that this indeed happens within our 
microscopic model and intrinsically complex supercon- 
ducting order parameter emerges naturally (of course, 
there is always an arbitrary overall phase). Such super- 
conducting state spontaneously breaks the time-reversal 
symmetry and minimizes the ground state energy for a 
range of coupling constants, which we determine below. 
For reader's benefit, we note here that an interesting and 
different possible time-reversal symmetry breaking order 
parameter, involving s- and <i-wave coexistence, was con- 
sidered in the context of pnictides in Ref. [26[ Finally, 
if one of the Josephson-like couplings between the bands 
is much smaller than the other two, one intuitively ex- 
pects that yet another form of the order parameter will 
appear: three gapped bands with a relative minus sign 
between the stronger-coupled ones. We show below that 
all of these possibilities are realized in different parts of 
the phase diagram of the microscopic model. 



II. THE MODEL AND ITS GAP EQUATIONS 

We start with a Hamiltonian which is a straightfor- 
ward generalization of the single-band BCS theory. Our 
model therefore bears all the birth-marks of the origi- 
nal - restriction to weak coupling, omission of many de- 
tails concerning band structure and dynamics of interac- 
tions, etc. - but shares some of its virtues as well: broad 
generality and simplicity which allows for analytic treat- 
ment. More realistic considerations would basically re- 
sult in various quantitatively important but conceptually 
straigthforward "decorations" of this simplified Hamilto- 
nian, which we now write down in its reduced form: 

+ E G™^c^c% + h.c., (1) 

ij,k,k' 

where the i and j are band indexes (they run from 1 

to 3) and for the moment we assume G% = (i.e. no 
intra-band pair-scattering). This simplifies the calcula- 
tions significantly and is justified by the following rea- 
soning. If we are to include the intra-band terms, there 
will be a finite critical strength for G% (k,k'), below 
which superconductivity cannot exist (for repulsive in- 
teractions). Above this threshold, when the supercon- 
ducting state is already present, the intra-band terms are 
irrelevant for the structure of the order parameter, which 
is entirely determined by the inter-band pair-scattering. 
This argument, however, has to be applied carefully (see 
below). The Josephson-like term G\ (k, k') is separa- 
ble and has the usual square-well form. We also assume 
identical parabolic two-dimensional (2D) bands. As we 
will see, in the gap equations it does not really matter 
whether we use hole or electron bands or some combina- 
tion. So our results apply for all of these cases, although, 
of course, the precise dynamics that produces supercon- 
ductivity in iron pnictides is most likely directly tied to 
its semimetallic character and the presence of both hole 
and electron bands at the Fermi level. This general na- 
ture of our results is the consequence of the simplified 
model and the relatively restricted set of question we are 
asking (for example, our focus is on the structure of the 
order parameter) . We can think of Eq. ([I) as describing 

effective low-energy theory and G% as phenomenolog- 
ical parameters in which we have stored all the details 
about the realistic band structure and the high-energy 
processes. 

Now, after we define mean field averages 

AL = - E ^(MO^Wk't) ( 2 ) 
and introduce Bogoliubov-transformed fermionic opera- 

(ii) 

tors, by using the properties of G 2 and following the 
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usual algebra, we get a set of three gap equation: 



,19,20 



( Wn ) 2 + {Eof 



> (3) 



where E % = -J (£ 1 ) 2 + (A 1 ) 2 , to n are the fermionic Mat- 
subara frequencies and ujc is high-energy cut-off. 

Despite the apparent similarity with the single-band 
BCS theory, these non-linear gap equations are consid- 
erably more involved, and do not allow for analytic so- 
lutions in the general case, even at T — 0. To achieve 
some progress, we need to simplify the model even fur- 

(23) 

ther. Let us start with two bands, coupled via G 2 , 
and then gradually turn on their couplings with a third 
band. We also assume these two new coupling constants 



( 1 2 ) [ i ' i 

G 2 and G y 2 °' to be equal. Thus, we reduce the three 
generally different couplings to two and introduce dimen- 
sionless constants: 

A< 12 > = N(0)G ( 2 12 \ A 13 = N(0)G { 2 13 \ X^ = N(0)G { 2 23) ; 

A (12) = A (13) = A) A (23)= r? . A ,77>0, 

where we have denoted the density of states (DOS) on 
the Fermi level as N(0) (identical bands!). With these 
simplifications, we are finally ready to make some ana- 
lytic progress and gain some insight of the physics of our 
model. 



-.(13) 



III. CRITICAL TEMPERATURE AND T 
REGION 



We now proceed by linearizing Eqs. ([3]) in the region 
T « T c and |A| < T. The problem then reduces to 
finding the eigenvectors and eigenvalues of a 3 x 3 matrix. 
The possible order parameters are proportional to the 
eigenvectors, and the eigenvalues determine T c . In this 
case the Eqs. © are equivalent to 



'0 A A\ /A 1 
/ | A rj A 2 I = 
,A v \ A 3 




(4) 



with I = 7 In (2uic/^T c ) > (7 is the Euler constant). 
Solving this matrix equation gives us three real eigenval- 
ues 

Si = -I V , I -(v-V8X^ + ^), ^ +v /8A 2 +f7 2 ),(5) 
and their corresponding eigenvectors 




For A, 77 > there are two negative eigenvalues and 
accordingly two possible order parameters. If we fix 77 
and gradually increase A from zero it is easy to see that 
eigenvalues 5i and S 2 cross at the point 77 = A. The Ai 



is obviously the two-gap solution and has higher T c for 
A < 77. The other possibility is a three-gap superconduc- 
tor with a relative minus sign between those bands that 
experience stronger coupling. Thus, the order parameter 
can be chosen to be real along the entire T c line. Pre- 
cisely at the crossing point, the eigenvalues are degener- 
ate and there any superposition of the two eigenvectors 
is also a legitimate order parameter. This degeneracy 
is a consequence of linearizing Eqs. Q and it leads to 
the possibility of complex A, with non-trivial phase dif- 
ference between the components. One example is the 
A cx {1, e~^ , e~~^} - the Ginzburg-Landau theory of 
this particular state was constructed and studied in Ref. 

Once T is below T c and one enters the superconducting 
state, we expect the complex order parameter to emerge 
as a competitor to the real one within a finite region, 
as opposed to a point at T = T c . Because of the 2 •<-> 3 
symmetry in the gap equations, we will look for solutions 
that satisfy the condition |A 2 | = |A 3 |. We can write the 
ordinary two- and three-gap order parameter as 

Ai = {0, -1,1)3, A 2 = { -0,1, 1}A (7) 

and introduce intrinsically complex, time-reversal sym- 
metry breaking order parameter of the form 



{-k, e 



lip 



-Z(p 



}fi. 



(8) 



In Eqs. and Eq. (JS| 3, 9, A, k, ft and ip are all real 
variables that parametrize the order parameters on the 
different bands and are to be determined self-consistently 
by solving the gap equations. 

To proceed analytically we expand the Eqs. ((3|) for 
(T c — T) /T c to second order in the magnitude of the order 
parameter. The two-gap solution obviously follows the 
BCS behavior: 



1 , ,2w c S 2 
- =7ln( — )-Po^2- 
V T 2 



(9) 



where we have simplified the notation by introducing a 
new constant f3 — 7((3)/8ir 2 . The three-gap order pa- 
rameters have more complicated behavior: 



«) h <^)4-g-i=o, 

, ( 2oj c A 2 
_= 7 l n (_)- A) 



(10) 



and 



,2w, 



Xk 3 



A K (l-. 2 )ln(^) + — -^ = 



CO sip 



2A ! 



ttT 
1 



V 



X 



1 ( 2UJc \ R ^ 

- =7ln( — )-Pq— 2 -. 



(11) 



We now solve these equations numerically and obtain 
Ai, A2 and A3. The time-reversal symmetry breaking 
(TRSB) solution exist only in the narrow interval A € 
(A c i, A C 2), where A c i, X c2 — > r\ for T — > T c . For A smaller 
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than A c i, the complex order parameter reduces to the 
two-gap one, while for the coupling parameter bigger 
than A C 2 only the trivial solution exist for A3 (see Fig. 
©)• 



1.00 



1.03 




FIG. 1. (Color online) Plot of re (solid line) and cos ip (dashed 
line) for A3 at T = 0.95T C . This OP exists as a distinct 
solution only at a vicinity of the degenerate point A = 77. On 
the left it merges with Ai and on the right crosses A2 at 
cosip = 1 and disappears. The interval is asymmetric with 
respect to the degenerate point. 

Using Eqs. (|9|), (fTO]) and (fTTj) we can proceed to cal- 
culate all the thermodynamic quantities of interest. We 
follow a simple route - solve Eqs. (|T0|) for 6(T,T C ) and 
Eqs. (HI|) for k(T, T c ), ip(T, T c ) (their structure allows it). 
After that we construct three single-variable Ginzburg- 
Landau theories (assuming uniform solutions), minimize 
the free energies and compare the results: 

•F = c^| 2 + §|V#, ^n = -^, (12) 

where \ip\\ = S, j^l = A and = ^- To do that we 
write the interaction part of the Hamiltonian as 



(«)-(»*) t-»t M 
2 L kt -kj." 



h.c. 



(13) 



and in Eq. (|13[) we have introduced auxiliary mean field 
averages d % = — X)k( c -kl c kt)- U sm g the definitions of 
{A 1 } we can write equations for d l : 



A 1 = 



A 3 = 



A (d^d 2 ), A 2 = -i-(Ad 1 +7 7 d 3 ), 



N(0) 
1 



N(Q) 



{\d 1 +i 1 d 2 ), 



N(0) 

which can be solved for {cf }: 

d^N(0)-^ 1+ ^ A2 + ^ 
d 2 = N(0) 
d 3 = N(0) 



2X 2 




77A 1 + A(-A 2 


+ A 3 ) 


2A?? 




-■qA 1 + A(A 2 


- A 3 ) 



Now we can obtain connected-diagrams expansion for the 
free energy in orders of d % . The second and the forth- 
order terms come from the expressions 

,1/T ,1/T 

Ti oc / dn / dT 2 (T T J^ mt (T 1 )J^ mt (T 2 )), 
Jo Jo 

,1/T ,1/T 
J4 CX / dTl.. dT4,(T T Jf int (Ti)..J#i nt (T4)). 
JO JO 



Using Eq. ([13)) for J^„ t , properties of the electron 
Greens' functions Sf 1 -? ~ 5y and expressing the d l, s via 
the A J 's we eventually get F in the form of Eq. (|12p . 
We have to calculate the pre-factors so we can compare 
the different order parameters. Expressions for a± (T, T c ), 
a 2 (T,T c ,9)), a 3 {T,T c , K , V ), 0i(T,T e ), fo{T,T c ,6)) and 
^3(T,T C , K,(p) are straightforward but tedious to obtain 
and very unwieldy, so we just report the results for J 7 . Far 
from the symmetric point A = 77 the solution with highest 
T c remains stable. On Fig. @ we show the comparison 
between different J-'s for T — 0.95T C at the vicinity of 
the symmetric point - the different solutions which are 
degenerate at T c split, and the complex order parame- 
ter has lowest free energy. This, however, remains true 
only in a relatively small interval around the line A = 77, 
which as T — > T c reduces to a point. This interval is 
asymmetric and considerably smaller on the A < 77 side. 
The different phases appear to be divided by a continu- 
ous phase transition on the left and a first order phase 
transition on the right side. 







1.00 



1.03 
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A/77 



2\rj 



FIG. 2. (Color online) Comparison between the Tsc — J~n 
for Ai (green), A2 (blue) and A3 (red). We show calculation 
on the right side of the A = 77 point for T = 0.95T C . The 
complex OP minimizes the free energy in a small interval. 
The transition between A3 and A2 is discontinuous. 



IV. LOW-TEMPERATURE REGION 

Now we will concentrate on T — properties of the 
model. To distinguish the parameters from the finite- 
temperature case we use subscript "0". To find S , 9o, 
A , k and f2 we have to solve T = version of Eq. 

The two-gap solution leads to identical gap magnitudes 
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which obey the BCS relation: 

S a 2u Jc e~ 1 '\ 



(14) 



it is 



Since A does not enter the gap equation for Ai, 
always a solution, irrespective of the ration rj/X. 

For the real three-gap solution A2 there are two un- 
knowns to determine. We get the following equations for 
them: 



1> 



-(7 - 26»o In O - 2 = 0, 



A = 2u c e 



-0o/2\ 



(15) 



Eqs. (fTSj) can be solved numerically and there is always a 
non-zero solution for 9q. That means that A2 is solution 
for all (non-zero) values of A. Looking at Ao, however, we 
see that it is strongly suppressed for A — ¥ 0, which is to 
be expected since at A = the only non-trivial solution 
is Ai. 

For the complex three-gap order parameter we get 
equations 



Ak ln(«o) - ko 

COS (fio 



X 2 



Xq 



0, 



n = 2 Wo e- 1 /i. 



(16) 



These equations have non-trivial solutions for A € (0, A c o) 
where A c o > r\ (see Fig. ((3J). At the point A = 77 the 
order parameter has the completely symmetric form 



a symm 
^3 



{1, e 



: }fio, 



(17) 



which is easy to understand if we consider the gap equa- 
tions. At this point So = f2o- 




FIG. 3. (Color online) Plot of no (solid line) and cost^o 
(dashed line) for A 3 at T = 0. This OP exists only for A e 
(0, A c0 ) where A c o > V- 

Now we want to see which order parameter is the ac- 
tual ground state for different A. We calculate the differ- 
ence between the superconducting and the normal state 
energies for the different A: 



£sc - £ 



N 



/zAy* A > - ^fs\^ - (iN p \y FS ) 



The normal and superconducting state kinetic energies 
are respectively: 



KB 



N 



ICS 



A 



£ 2& 



(18) 



Ml! 



Converting the sum into an integral gives 



K,£ K - K,£ N = 2 J2 N (°) 
N(0)J2 ' -A' l" sinh 





1 ( fc_ 



m 2 



.(19) 



The mean field average of the potential energy in the 
normal state is zero, and for calculation in the supercon- 
ducting state we again use d % . Then the potential energy 
can be written as 

VE K . ~ V£ N = Xd x *d 2 + Xd^d 3 + rjd 2 *d 3 + fc.c.(20) 

Using the expressions for d % in the potential energy 
formula we get 

-77IA 1 ! 2 + A 2 (-|A 2 | 2 + |A 3 | 2 ) 



vs A -vs N = N(oy 



A 2 r? 



•JV(0) 



Ar/(A 1 (A 2 * + A 3 *) + A**(A 2 + A 3 )) 
X 2 rj 



Combining Eqs. (JT5]) and (|2U1) we can compute the 
energies for the different possible ground states and com- 
pare them. Let us start with Ai - since So follows the 
BCS behavior we get the standard result, multiplied by 
2 (two bands): 



£ 



N 



-iV(0)S, 



—2 

o- 



(21) 



Similar calculations for A2 and A3 give the energy dif- 
ference as a function of 9q and Aq, or kq and Hq: 



£ A 3 



£ 



N 



-N(0)(l + ^jA 2 , (22) 



(23) 



K 



^(0) l-- + -^(l--+21n«o) )% 



V 



N(0) 



AXkq cos ipo 2 sin <po 



2X 2 



V 



n 2 

S2 . 



Now we are able to compare the different solutions - 
the result is shown on Fig. (|4|). On the left side of r\ 
we see that the time-reversal symmetry breaking order 
parameter is the ground state. It converges from below 
to the two-gap solution as A, #0 — * 0. At the symmetric 
point A = r\ we can use the BCS result for both solutions 
and £^ a — £m — 3/2(£^ — £n) (three vs two gaps). On 
the right side there is a Quantum Phase Transition at 
some X cr > 77, where the complex and the real three-gap 
states merge (cost/? = 1). Beyond this point A3 ceases to 
exist. 
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0.9 1 1.1 

FIG. 4. (Color online) Comparison between the £sc — £n 
for Ai (green), A2 (blue) and A3 (red). The first one is never 
a ground state for A 7^ 0. The energies for A2 and A3 merge 
for some X cr > r\. 
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FIG. 5. (Color online) Suggested phase diagram of the three- 
band model. There are three possible superconducting (SC) 
order parameters (OPs). The line separating the TRSB and 
real three-gap OP is most likely first order phase transition 
line. 



V. PHASE DIAGRAM 

On the basis of the above results we suggest that our 
model has the phase diagram depicted in Fig. (|5|) . There 
are three superconducting order parameters, stable in dif- 
ferent regions, separated by two critical lines. On the left, 
different superconducting states appear to be separated 
by continuous transition and on the right by first order 
one at finite temperature, and continuous one at T = 0. 
There is a possibility of observing two different supercon- 
ducting states in a single system and the transition be- 
tween them, tuned by the temperature. In the context of 
this phase diagram we consider the case of iron pnictides 
- bands 2 and 3 can be thought as the hole bands at the 
r point, which are strongly coupled to one of the electron 
bands at the M = (ir, tt) points. If the renormalization 
group arguments apply, G^ h ^ are enhanced and G^ 1 ^ 2 
is suppressed by the same high-energy electron-hole pro- 
cesses. Then the appropriate regime for the pnictides is 
A > 77, on the right side of our diagram. The existence of 
time-reversal symmetry breaking order parameter is not 
excluded, but is unlikely for the optimally doped com- 
pounds, given the relative narrowness in which it is stable 
for A > rj. However it may be present in the overdoped 
materials, for which the inter-band interactions are sup- 
pressed, due to the significant deviations from perfect 
nesting. This complex order parameter, in general, en- 
tails the existence of local magnetic fields at edges and 
around impurities, and likely domain structure^,. These 
effects may provide the best way of observing such state 
and its broken symmetry. 

Before we proceed, let us comment on two obvious de- 
ficiencies in our model, which seemingly prevent us from 
applying the results we have derived thus far to the iron 
pnictides. First, we have completely neglected the intra- 
band pairing terms. Once these terms are included the 
calculations become considerably more involved and it 
is difficult to proceed short of pure numerics. However, 
we believe that our phase diagram is qualitatively cor- 



rect even in that case, since, as already explained, these 
terms do not play a role in determining the structure of 
the order parameter (provided that the superconductiv- 
ity is still possible) and only change the numerical val- 
ues of various results (T c , for example). The validity 
of this argument is limited, however, and the intra-band 
terms have an important role to play in the case of several 
competing channels (s and d— wave, for example) which 
are affected differently by these terms. There are several 
studies for pnictides suggesting such competitio n 1 1 i 29 i 30 . 
But as long as the most isotropic superconducivity re- 
mains the leading instability, it will be realized without 
any mixing from the sub-leading channels (for s and d 
mixing see Ref. OHI) and our results apply. In case the 
system is driven to a nodal state by the intra-band re- 
pulsion, the frustration due to the inter-band terms can 
again lead to a development of complex order parame- 
ter, but we leave this question for further studies. Sec- 
ond, we have restricted ourselves to a three-band model, 
whereas in pnictides there are generally four active bands 
participating in the superconductivity (see, for example, 
Ref. \m. It is a valid question if adding another band 
will completely suppress the complex order parameter. 
To address it let us remind the reader the tight-binding 
calculation 16 , which indicates that, for the case of pnic- 
tides, the pair-scattering terms between the second elec- 
tron band (whose existence we have neglected) and the 
hole bands are at least an order of magnitude smaller that 
Qeihi ( ana logous to our G^ 3 and G| 3 terms). Coupling 
between the electron bands, however, generically will be 
of the same order as G hlfl2 (our G^ 2 )- This means, in 
practice, that the phase of the gap on the second elec- 
tron band will (almost) entirely depend on the gap on 
the other electron band. Thus the relative phase between 
the gaps opening on the hole and the electron bands will 
still be determined by the three-band calculation. These 
(somewhat naive) arguments allow us the hope that our 
model, despite its simplicity and numerous assumptions, 
is relevant for the iron pnictides. 
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One condition of particular relevance to the pnictides 
is the condition for existence of superconductivity itself 
- in the two-band model it is Gf > U, (G^ 1 = Gf = U 
is the intra-band pairing). For the real three-band order 
parameter in the limit A 3> r/ this condition becomes 
i.e. it is somewhat relaxed. 



two-gap solutions (previously A solely was driving the su- 
perconductivity on the third band). The equation for T c 
becomes: 



(26) 








A 2 ) = - 


[ A 2 


A 3 / 


\A 3 



VI. JOSEPHSON-COUPLED TWO-GAP S ' 
STATE AND SINGLE-GAP S STATE 

We can use the model and the results derived so far to 
study a different problem - a two-gap s' state coupled via 
Josephson junction to an ordinary s-wave superconduc- 
tor. This is a situation of real experimental relevance, in 
light of the recent experiments demonstrating Josephson 
effect between Pb and an iron-pnictide superconductor 32 ; 
a theoretical background is explored in^ 3 - - — . The tun- 
neling of Cooper pairs in this case would like to align 
the phases of the two (distinct) superconductors - we 
can model this by introducing negative coupling constant 
A — > —A (and |A| <IC 77 - weak coupling). It does not take 
one much time to realize that the equations for this model 
can be made identical to the ones for the previous model 
by a single sign flip in the three-gap order parameters. 
For example, the real solution A2 now becomes: 



{ -9 ,1, 1}A4{ 6 ,1, 1}A. 



(24) 



After this change the phase diagram is identical to that 
on Fig. (|5|). The model is still frustrated, but the frus- 
tration is resolved in a different manner - the non-trivial 
phase angle now brings the 77-coupled gaps closer (instead 
of further away) to the third gap. This is easy to see at 
the completely degenerate point A = 77 (Fig. ((6])): 



X symm 
^3 



{1, e^ , e-*}n. (25) 
The model with negative Josephson junctions still does 




FIG. 6. Schematic representation of the TRSB order param- 
eter in the case of three positive inter-band couplings (left) 
and two negative and a positive inter-band couplings (right). 
The frustration is resolved in a different but related way. 

not give us two independent superconductors. To achieve 
this we add intra-band attraction on the weakly-coupled 
band. Now even at A = we have two different su- 
perconducting states - a single-gap s-wave and a s'-wave 



where we have parametrized the intra-band attraction 
as a fraction b of 77. For the experimental set-up of a 
conventional low-temperature superconductor coupled to 
iron pnictide sample we expect b < 1. The eigenvalues 
and eigenvectors now are: 



Si 



-fa(l-6)T 0U 2 +r; 2 (I + 6) 2 ), (27) 




??(l+fe)±^/8A 2 +T) 2 (l+&)~ 

f 
f 



(28) 



Again there are two possible order parameters. The T ci 
curves (Fig. ([7])) still cross but their crossing point is no 
longer at A = 77. It moves to the left, which is easy to 
understand - the T c line for the three-gap OP goes to a 
finite limit rather than zero for A — > (single-gap SC, 
courtesy of the non-zero b). 







0.25 0.5 0.75 



L 



A/17 



FIG. 7. (Color online) Comparison of the T c for Ai (blue) 
and A 2 (red) for 6 = (solid line), b — 0.5 (dashed line) 
and 6=1 (dotted line), n is fixed. For small A the two-gap 
solution is the first to appear for 6 < 1. At the crossing point 
there is degeneracy and complex A is possible. 

We again expand the gap equations in the vicinity of 
T c up to second order in the magnitudes of the supercon- 
ducting gaps. The two-gap solution is unchanged, and 
the real three-gap solution equations become: 



(Afl-,)-(Afl 3 -,) 2A + ^ 3 W**- 



(V ~ A# 3 ) 



2A + bnO 3 



-1 = 0, 



7 (2A + brj6) ln( 



~ttT 



(2A + b V d 3 )l3 



A 2 
T 2 ' 



(29) 
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The complex order parameter equations are: 



Or 



2u r 



Ak 3 



A/c(l - «f ) ln(— 7) + 



t]Ak 



costp 



ttT 
tjXk 



r] (br) 2 + A 2 ) 



0. 



1 2w ft 2 

2(6^ + A 2 )' ^ =7ln( ^T ) ~ /3o 7f (30) 



For 6 — > these equations reduce correctly to the 
inter-band-couplings-only case. We derive single- variable 
Ginzburg-Landau free energy, and then minimize it with 
respect to A and tt. The comparison between the differ- 
ent solutions for b = 0.5 is shown on Fig. (JHJ. The result 
is very similar to the b — case but the region in which 
the complex order parameter dominates is smaller. The 
reason is that the real three-gap solution's free energy is 
pushed down by the intra-band term. However, for b — > 1 
the region again expands as the crossing point is pushed 
closer to A = 0. For b > 1 the real three-gap solution 
minimizes the J- for all A, at least for T « T c . 




0.71 



0.73 



A/77 



FIG. 8. (Color online) Comparison between the Tsc — J~n 
for Ai (green), A2 (blue) and A3 (red) for b = 0.5. We show 
calculation for T = 0.95T C at the vicinity of the T C i crossing 
point A ~ O.7I77. The interval for which the complex OP 
minimizes the free energy is smaller but still exist. 

Now we discuss the T = line of the phase diagram. 
The gap equations for A2 and A3 are now: 



9o-j0 0- - +r?ln0o)) - 20 o ln0 o -2 = 0, 



A = 2uj c e 



-9 a (l+b V In e )/(2\+br,e Q ) . 



and 



Ak ln(/c ) - k 



r)K 



b-rK (l - r]\nn ) = 0, 



cos — ^—^-(1 — b + brj In k ), Qq = 2u>ce 1 ^ n . 
2 A 

We compare the energies for the different possible ground 
states on Fig. Consistent with the results from the 
T c region, the complex solution is still the ground state 
on the left of some A cr . However, because the energy of 
A 2 is pulled down for non-zero b, the transition a weakly 
first-order one and A cr is on the left of the T C i crossing 
point. 



I 

U 
00 

f*3 



0.5 



0.7 
A/77 



0.9 



FIG. 9. (Color online) Comparison between the Esc — £n for 
the different order parameters. Here b = 0.5. The energies 
for A2 and A3 cross for A cr w 0.65. 



With increase of b, A cr moves to the left, but there is 
always a region (confined to lower and lower tempera- 
tures and smaller and smaller A as b goes up) in which 
the complex solution is the preferred order parameter. 



VII. CONCLUSIONS 

In summary, we have considered a simple microscopic 
model, with three bands coupled via repulsive pair- 
scattering interactions, which is relevant for the recently 
discovered iron-based family of high-temperature super- 
conductors. We have constructed the phase diagram of 
this model and discussed its overall features. Generally, 
we find three possible superconducting order parameters, 
one of which breaks the time-reversal symmetry in order 
to relax some of the frustration intrinsic to the three (or 
odd) band case. The conditions for such exotic state are 
rather strict and it seems unlikely that this order param- 
eter would be observed in the optimally doped iron pnic- 
tides. However, this state may be realistically present 
in overdoped samples, if the doping is carefully tuned to 
the range of optimized frustration. While quantitative 
aspects of our results are bound to be sensitive to the de- 
tails of the band-structure and the accompanying orbital 
character of each individual iron-pnictide material - the 
details which are not part of our model - the overall qual- 
itative features reported in this paper are expected to re- 
main relatively universal. Experimental observation of a 
time-reversal symmetry breaking superconducting state 
is perhaps the best we can hope for in linking an s' super- 
conductor to some broken symmetry and would represent 
arguably the strongest confirmation yet of the basic pic- 
ture which places the repulsive, purely electron-electron 
interband interactions at the heart of iron-based high- 
temperature superconductivity. Furthermore, we have 
also considered the case of Josephson-coupled two-band 
s' SC and a single-gap s SC. Again, there is possible time- 
reversal symmetry breaking state, although frustration in 
that case is relieved in a different (but related) manner. 
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